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The giant monopole, dipole and quadrupole responses in 40Ca, 90Zr, ^^Osn 

' and 208p)-, g^pg investigated using linear response treatment based on a 

' stochastic one-body transport theory. Effects of the coupling to low-lying 
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surface modes (coherent mechanism) and the incoherent mechanism due 
to nucleon-nucleon collisions are included beyond the usual mean-field de- 
scription. We emphasize the importance of both mechanism in the frag- 
mentation and damping of giant resonance. Calculated spectra are com- 
pared with experiment in terms of percentage of Energy- Weighted Sum- 

' Rules in various energy regions. We obtained reasonable agreement in all 

\ 

cases. A special attention as been given to the fragmentation of the Gi- 
ant Quadrupole Resonance in calcium and lead. In particular, the equal 
splitting of the 2+ in '''^Ca is correctly reproduced. In addition, the ap- 
pearance of fine structure in the response 208p|j jg partly described by 
the calculations in which the coherent mechanism play an important role. 
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I. INTRODUCTION 



Recent development of high resolution experiments offers the possibility for a deeper un- 
derstanding of collective motion in quantum fermionic system like nuclei. These experiments 
enable to determine the fragmentation of the nuclear response with a very high resolution up 
to few keV Understanding of the fine structure in the nuclear collective response, its 

fragmentation and damping mechanisms constitute a major challenge for theoretical models 
1^-^. One possible avenue is the development of quantum transport theories for nuclear 
dynamics |^,|TU . 



In dynamics of nuclear motion, one usually distinguishes damping due to the coupling 
to the external and the internal degrees of freedom. The former one gives rise to cooling 
of the system by evaporation of particles, while the latter one leads to the dispersion of 
the well ordered motion through mixing with the internal degrees of freedom. In the latter 
case, one can again distinguish (i) the Landau damping due to spreading of the collective 
modes over non- collective particle-hole (p-h) excitation, (ii) the coherent mechanism due to 
coupling with low- lying surface modes and (iii) the damping due to coupling with the 

incoherent 2p-2h states usually referred to as the collisional damping |[12|jl3| . 

Most investigations of the nuclear response carried out so far are based on either the 
coherent damping mechanism or the collisional damping. The coherent mechanism is, par- 
ticularly, important at low temperature, and accounts for the main feature of fragmentation 
of the response |p|,p!li[l4HT7|| . On the other hand, the collisional damping is relatively weak 



at low temperature [|l8l, but its magnitude becomes large with increasing temperature, as 



shown in recent calculations ||T9H22|. In this work, we carry out investigations of nuclear 
collective response on the basis of a one-body stochastic transport theory, which incorpo- 
rates both the coherent mechanism and the collisional damping in an consistent manner 



as demonstrated in p3| , p4| (also see p5| , |26[| ). We calculate the giant monopole, dipole and 
quadrupole responses in '^'^Ca, ^°Zr, ^^''Sn and ^°^Pb , and compare the results with experi- 
ment in terms of Energy- Weighted Sum- Rules distribution. We find that both mechanisms 
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play important role for a proper description of the fragmentation and the damping of giant 
resonance excitations. 

In section II, we present a brief description of the linear response treatment of collective 
vibrations based the stochastic one-body transport theory . In section III, we discuss the 
details of the calculations and present the results and comparison with data in section IV. 
Finally, we give the conclusions in section V. 



II. LINEAR RESPONSE BASED ON STOCHASTIC TRANSPORT THEORY 

1. Stochastic transport equation 
In the stochastic transport theory, temporal evolution of the fluctuating single-particle 



density matrix p{t) is determined by |10 



zh^^pit) - [h{p),p{t)] = Kj{p) + 6K{t) (1) 

where the left hand side corresponds to the mean-field evolution in terms of the self-consistent 
mean-field Hamiltonian h{p) expressed in terms of the fluctuating density, and the right hand 
side arises from the correlations due to residual interactions. The first term K/(p), which 
is usually referred to as the binary collision term, describes the coupling of single-particle 
excitations with more complicated two-particle two-hole states. It can be expressed as, 



Klip) = f k t/i2(t, s)F,2{s)uUt, s)] ds (2) 



where Ui2(t,s) represents a product of the mean- field propagator Uu = Ui ® U2 with 
U{t,s) = exp (^-i/h J\{p{t'))dt'^ and 

F12 = (1 - Pi)(l - P2)v/hP2 - PiP2v{l - Pi)(l - P2)- (3) 

Here pip2 Represents the anti-symmetric product of the single-particle density matrices and 
V denotes the residual interactions. As seen from eq.(2), the collision term, in general, 
involves memory effects due to the time integration over the past history from an initial 



time to to the present time t. The second term in the right hand side of eq.(l) is the initial 
correlation term, 

6K{p)=Tr2[v,6au{t)]. (4) 

where Sai2{t) = Ui2(t,tQ)S(Ti2{to)Ui2(t,to) represents the propagation of the initial correla- 
tions from to to t. In the stochastic transport description, the initial correlations Scri2{t) are 
treated as a Gaussian random quantity. Consequently, the initial correlation term 6K{t) 
has a Gaussian distribution characterized by a zero mean and a second moment, which can 
be determined in accordance with the fluctuation-dissipation relation of the non-equilibrium 
statistical mechanics. 

In the stochastic transport description, dynamical evolution is characterized by con- 
structing an ensemble of solutions of the stochastic transport eq.(l). In this manner, the 
theory provides a basis for describing the average evolution, as well as, dynamics of density 
fluctuations. In the semi-classical approximation, a number of applications have been carried 
out for description of multi-fragmentation in heavy- ion collisions p7| , |28| . Furthermore, as 
demonstrated in recent publications I^S],^, the stochastic evolution involves a coherent dis- 
sipation mechanism arising from the coupling of single-particle motion with the mean-field 
fluctuations. This can be shown by considering the average evolution of the single-particle 
density matrix p{t) = 'p{t). The ensemble average of eq.(l) is calculated by expressing the 
mean-field and the density matrix as h{p) = h{j)) + 5h{t) and p{t) = p{t) + Sp{t), where 
6h{t) and Sp(t) represent the fluctuating parts of the mean-field and the density matrix, 
respectively. Then, the ensemble averaging yields a transport equation for the evolution of 
the average density matrix, 

^^Pit) - Hp),m] = Kjip) + Kcip) (5) 

where Kj{p) represents the incoherent collision term and the additional term arises from the 
correlations of the mean-field fluctuations and the density fluctuations, 

Kc{p) = [6h{t),6p{t)] (6) 
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and it is referred to as the coherent colhsion term. For small fluctuations around the aver- 
age evolution, the density fluctuations can be expressed in terms of RPA phonons, and the 
coherent term takes the form of a particle-phonon collision term. As a result, eq.(5) pro- 
vides an extended mean-field description, which goes beyond the extended Time-Dependent 
Hartree-Fock theory by including a coherent collision term into the equation of motion in 
addition to the incoherent one Q. 

2. Linear response based of extended mean-field theory 

In this section, we consider the small amplitude limit the transport eq.(5) and give a brief 
description of the linear response formalism including both the incoherent and the coherent 
damping terms. A details description of the formalism can be found in recent publications 

|2T|,g. 

The linear response of the system to an external perturbation can be described by consid- 
ering the small amplitude limit of the transport eq.(5). The small deviations of the density 
matrix 6p{t) = p{t) — po around a finite temperature equilibrium state po are determined by 
the linearized form of the transport eq.(5), 

in^5p - [/lo, 5p] - W + A, Po] = 5K\p) + 5K^(p) (7) 

In this expression A{t) = Aexp{—iujt) + h.c. is a one-body harmonic excitation operator 
containing a small imaginary part uj = uj + irj, and 6Kj{p) and 6Kc{p) represents the 
linearized form of the non-Markovian incoherent and coherent collision terms, respectively. 

The steady state solution of eq. (|^) can be obtained by using a development in terms of 
the RPA functions, 

exp(— icDt) + h.c. (8) 



In the following, we denote the average one-body density p by p. 



Sp{t)= Q+,po 



5 



where = J2x>o 0\ — z^^Ox. In this expression, 0\ and 0\ are the creation and 
annihilation operators of the RPA state A of energy hux , which are determined by the finite 



temperature RPA equations [2£ 



hujxOi=[ho,Q{] + hl 



(9) 



where h\ = {dh/dp) ■ p\. Substituting the expression (8) into eq.(7) gives rise to a set of 
coupled equations for the amplitudes z'^ and z^ coefficients that can be recast into a matrix 
form pT|, 



{huj — S {uj)) 



(10) 



where z^ and z~ are the amplitude vectors with components z^ and z^ , A is the forcing 
vector with components Ax = Tr[OA, A]pq and S (a)) denotes the self-energy matrix. In the 
small amplitude limit, the self-energy can be separated into the incoherent part and the 
coherent part, S (cD) = S/ (u) + Sc {(^)- According to the expression of the incoherent 
part, which also contains the RPA energy, is given by 



Am M 



LJxSx, + K^; 



UJ] 



-K 



A/1 



-uo 



(11) 



In the Hartree-Fock representation, the elements of the incoherent self-energy are given by, 

(12) 
and 



Kt; (Co) 



1^ {M\ [Ox,v] \^J) {tj\[Olv]\ki) ^^ 

'7 ^ A ■'^ijki 

4 tj^e ^ - ^^ijki 



Kt- {^) 



ly {H\ [Ox,v] \i3) {i3\\0,,v]\M) ^^^^^ 



(13) 



with MijM = (1 - ni){l - nj)nkne - ^^^^(l - nfc)(l - rig) and AeijM = Si + Ej - Sk - ei, 
where Si and rii denote energies and Fermi-Dirac occupation numbers of the single-particle 
states. The collisional self-energy is non diagonal, and therefore it introduces a coupling 
between different RPA modes through their decay channels, so-called collisional coupling. 
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In the following, we will neglect the non-diagonal part, which in general introduces a small 
correction to the strength distributions. 



According to the expression of the coherent self-energy is given by 

|2 



^ I <^\[Q,^^h{]\j > I 

— — —^—^>^^ij 



Xij 



huj — huJx — tj + Cj 



(14) 



^ I <t\[Qf,,hx]\3 > P 



Xij 



hCo -\- hu\ — ej + e 



where 



Mx,, = (iV° + 1)(1 - n;)n° - NXil - n°) 



(15) 



and denotes the finite temperature boson occupation factors for the RPA modes = 
l/[exp{huj\/T) — 1]. In general, the coherent self-energy is also non-diagonal, and it couples 
different RPA modes. Here, we neglect this coupling and retain only the diagonal part. 
The coherent mechanism, which arises from coupling of the single-particle excitations with 
the mean-field fiuctuations in the stochastic transport theory, correspond to the coherent 
mechanism described in and its finite temperature extension using the Matsubara 



formalism in 14.15 



We can deduce the response function associated with an excitation operator A , by 
calculating the expectation value < A >= Ti A6p(t) with the help of the expression (8). 
The corresponding strength distribution is obtained by the imaginary part of the response 
function and it can be expressed as. 



1 



TT 



S{u) = —Imi A*, A ){hCo-J:{Co)) 



,-1 



/ 



A 

-A* 



(16) 



The strength function includes both damping mechanisms, i.e. the collisional damping due 
to coupling with the incoherent 2p-2h states and the coherent mechanism due to coupling 
with a low-lying phonon and p-h states. 

In our previous studies, we investigated the nuclear collective response in the basis of 
the incoherent damping mechanism. We found that at low temperature, in particular for 



light and medium weight nuclei, the incoherent damping mechanism has a sizeable influence 
on the strength functions, and it becomes more important at higher temperatures. On the 
other hand, in particular for heavy nuclei, the coherent mechanism due to coupling of giant 
resonance with phonons plus p-h states, plays a dominant role for describing the properties 
of cold giant resonance. In this paper, we want to clarify the relative importance of the 
incoherent and the coherent mechanisms in collective response in cold nuclei. For this pur- 
pose, we present three different calculations by incorporating only the coherent mechanism, 
only the incoherent mechanism and including both mechanisms in to the calculations, and 
compare the results with the experimental data. 



III. DETAILS OF THE CALCULATION 



A. RPA calculation 



In order to obtain the solution of equation (|T0]), we first solve the RPA equation @ in 
a discrete basis. In order to account partially for the states in the continuum, particle and 
hole states are obtained by diagonalizing the Hartree-Fock Hamiltonian in a large harmonic 
oscillator representation which includes respectively 12 major shell for ^°Ca and 15 major 
shells for other nuclei. We use a fixed imaginary part for the forcing frequency 1] = rj^ = 0.5 
MeV. The Hartree-Fock and RPA calculations are performed using the effective Skyrme force 
SLy4 We use the standard excitation operators for isoscalar and isovector resonances 
for L > 1 (for a review see 



(17) 



^LM — -J l^i=i r i Ylm - ^ 2^i=i r ^ Ylm 
where Ylm are the spherical harmonics, and for isoscalar giant monopole resonance, we 
employ 

Aoo = ^j:rlY,o (18) 



The energy-weighted sum rule (EWSR) is given by, 

m, = Y.huj,\{0\A\\)f (19) 

A 

When the states A are specified in the RPA, it can be calculated from the Hartree-Fock 
ground state according to mi = | < [F\ [H,F]] >o. For the Skyrme interactions, it leads 
to the following expression, 

^GMR _ 27i2 z2 / 2\ 



(20) 



^GQR _ 50 Z2 



47r 2m A ) HF 

where {j''^)hp denotes the root-mean square radius (rms) obtained from the Hartree-Fock 
ground state. In the case of the Giant Dipole Resonance, the Thomas-Reich Kuhn (TKR) 
sum rule is violated due to the non-local term in Skyrme forces, and the modification factor 
is giving by p2| , |33| . 



2m 



1 ^ . In 

tl{l + -Xl) + t2{l + -X2) 



1 



(21) 



where pn and pp are the neutron and proton one-body density. 

In the following, we compare the result of calculations with the experimental EWSR by 
employing the standard expression and parameters of the sum-rule |3^J3^. In the standard 
approach, the rms radius is approximated using a Wood-Saxon shape for the one-body 
density, which leads to the following expression, 

~ 'na' 



ws 5 ^\ 3 



1 2> 



(22) 



where Rq correspond to the surface position and a is the diffuseness. Different parameters 
used in the calculations are reported in table |. In table we compare the sum rules 
obtained by the parameterization ( ^2]) and k = and those results obtained from the 
RPA calculations of (|T^). The smallness of the difference insures the quality of the RPA 
calculations. 
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TABLES 



Nucleus 



Rq 



V /ws 



Ihf 



K 



(fm) 



(fm) 



(fm) 



(fm) 



^OCa |6| 
90Zr m 



3.65 
4.90 
5.55 
6.67 



0.55 
0.515 
0.515 
0.545 



3.49 
4.25 
4.71 
5.55 



3.40 
4.26 
4.70 
5.55 



0.160 
0.177 
0.176 
0.180 



TABLE I. Density shape parameters used in the calculations of the sum-rules. Hartree-Fock 
rms obtained with the Sly4 force are very close to those obtained with the Wood-Saxon parame- 
terization. 



Nucleus 



0+ (MeV. fm" 



1- (MeV. fm^ 



2+ (MeV. fm" 



40 



Ca 



10103 (9290) 



148 (144) 



10050 (9226) 



90 



Zr 



26664 (26025) 



330 (325) 



26523 (22556) 



120 



Sn 



38295 (37026) 



433 (427) 



38092 (26477) 



208 



Pb 



82633 (84572) 



738 (748) 



82197 (76034) 



TABLE II. nil sum rule obtained from standard parameterization of the rms radii. Sum rule 
calculated from RPA are reported in parenthesis. 
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B. Computation of self-energies 

1. Coherent mechanism. 

In order to incorporate the effect of coupling to surface modes, we calculate the RPA 
response for multipolarities up to L = 5. The coherent self-energy given by eq. (14) is then 
calculated by coupling collective states to low- lying states that exhaust at least 1% fraction 
of the EWSR. Energies and EWSR of collective modes used in the calculation are reported 
in table We note that, in particular, energies of collective 3^ states are overestimated in 
the RPA calculations. The percentage of the EWSRs are normalized to the RPA sum-rules. 
In the calculations of the matrix elements in eq.(14), we employ the full SLy4 interaction. 



Nucleus r E (MeV) % EWSR 



0+ 


17.6 


10.4 


0+ 


19.1 


12.5 


0+ 


20.7 


18.1 


0+ 


22.0 


18.0 


0+ 


24.4 


15.5 


1- 


16.7 


19.5 


1- 


17.8 


15.0 


1- 


18.6 


22.8 


2+ 


17.1 


74.9 


3- 


5.3 


9.5 


3- 


7.3 


10.8 


5- 


5.2 


3.7 
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90Zr 0+ 17.8 35.1 

0+ 18.6 17.1 

1- 15.7 50.1 

1- 17.9 14.9 

2+ 5.4 5.4 

2+ 16.3 13.5 

2+ 16.4 54.4 

3- 3.7 4.9 

3- 7.9 19.8 

3- 9.4 8.9 



^20Sn 0+ 16.8 30.7 

0+ 18.4 42.4 

1- 14.1 20.1 

1- 14.5 13.7 

1- 16.9 15.7 

2+ 6.0 3.5 

2+ 16.5 23.0 

2+ 16.6 24.0 

3- 3.8 6.7 

3- 6.9 7.2 

3- 7.6 18.0 

5- 7.0 2.7 

5- 8.3 1.1 

5- 9.5 5.8 
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208pb 0+ 

0+ 
0+ 

0+ 

1- 
1- 
1- 

2+ 
2+ 
2+ 
3- 
3- 
3- 
3- 
4+ 
4+ 
4+ 
4+ 
5- 
5- 
5- 
5- 
5- 
5- 
5- 



14.4 

14.7 

15.0 

16.4 

12.8 

13.2 

13.7 

3.7 

6.1 

13.0 

4.9 

6.7 

7.2 

9.1 

6.6 

8.9 

9.3 

9.6 

4.8 

5.7 

5.9 

6.7 

7.7 

8.2 

9.3 



33.5 
18.8 
19.5 
14.3 
13.7 
21.2 
13.8 
1.4 
10.2 
67.5 
21.9 
2.0 
5.7 
1.3 
4.0 
3.1 
3.1 
3.1 
1.3 
2.3 
1.1 
1.1 
2.5 
1.7 
1.6 
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5- 9.5 1.9 

TABLE III. Collective modes and associated EWSR obtained from RPA that are included in 
the calculations. 
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2. Incoherent mechanism 



We have shown in previous calculation that Skyrme interaction are not adequate to 
compute the incoherent part of the self-energy due to the presence of high momentum 



component (also, see |jT3[)- As in previous application |[T9|, we use a modified Skyrme 
interaction which is obtained by introducing a Gaussian cut-off factor in the matrix elements 
of the Skyrme force 

/ /32|<q2>|\ 

< %2\v\kl >=< ij\vs\kl > ■ exp ( j-^ 1 (23) 

in this expression, j3 describes an effective range of the interaction and the quantity < 
> provides a measure for the relative momentum which is defined by the relation < 
ij\6{r)\kl >< >=< zj|q^(5(r) + 6{r)c{^\kl > with r = ri — r2 and q = (pi — P2)/2. A 
quantitative discussion of the influence of (3 can be found in p2|. In the following, we use 
/3 = 1.4 fm. The size of the HO basis used to expend the particle and hole (HF) states have 
been chosen large enough to ensure a convergency of the results with in few per cent. 



IV. RESULTS 

We carry out calculations of strength functions for the giant monopole (GMR), the giant 
dipole (GDR) and the giant quadrupole (GQR) response at zero temperature. The resulting 
strength distributions are presented for "^^Ca (figure |l]), ^°Zr (figure |^), ^^°Sn (figure ||) and 
^°^Pb (figure ^). The left panels of figures show the result of RPA calculations (dashed 
lines) and the calculations performed by including only the coherent self-energy (thin lines) 
and only the incoherent self-energy (thick lines). The result of calculations performed by 
including both self-energies are shown on the right panel of figures (thick lines). 
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FIGURES 
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10 15 20 25 10 15 20 25 
oj (MeV) u (MeV) 



FIG. 1. Calculated strength distributions for GMR (top), GDR (middle) and GQR (bottom) 
in ^"Ca. Left: strengths obtained in the RPA (dashed lines), with the coherent mechanism (thin 
lines) and the incoherent mechanism (thick lines). Right: comparison between the RPA (dashed 
lines) and the extended RPA (thick lines), which includes both the coherent and the incoherent 
damping mechanisms. 
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5 10 15 20 5 10 15 20 

oj (MeV) Gj (MeV) 

FIG. 4. Same as figure |l| for ^ospb. 



A. Interplay between incoherent and coherent mechanism 

In order to quantitatively discuss the effects of different contributions, it is useful to 
compute moments of the strength in a given energy interval, 

^i-[s^,„,s™a.] = / {hujy S{w)duj (24) 

From these moments, we can define various mean energies Ei = mi/mi-2. An estimation of 
the spreading of the strength is given by the width F = \Jm2/mQ — {mi/mof'. In figure |^, 
the difference AE = Ei — Ei between the mean energy obtained in different calculations 
and the mean-energy calculated in the RPA is plotted as a function of the mass number for 
the different multipolarities. In figure, calculations by including the coherent mechanism, 
the incoherent mechanism and the coherent plus incoherent are indicated by dashed lines, 
dashed-dotted lines and solid lines, respectively. In figure ^, a similar plot is presented for 
the deviation of the width from the RPA response, AF = F — y^^"^ . In these calculations. 
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moments of the strength functions are evaluated over the energy interval 0-40 MeV. Prom 
the results of calculations, we can draw the following conclusions: 



— □(!) 

0.5 

0.0 - 
-0.5 7 
-1.0 -_ 

1.0 7 




-1.0 h. A 

50 100 150 200 
MASS 



FIG. 5. Variation of the mean-energy AE = Ei— Ei calculated in the energy interval 0-40 
MeV for GMR (top), GDR (middle) and GQR (bottom) as a function of mass number. Calculations 
performed by including the coherent mechanism, the incoherent mechanism and both coherent and 
incoherent mechanisms are indicated by dashed lines, dash-dotted lines and solid lines, respectively. 
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JIPJ^ 

FIG. 6. Variation of the mean-width AE = Ei — Ei calculated in the energy interval 0-40 
MeV for GMR (top), GDR (middle) and GQR (bottom) as a function of mass number. Calculations 
performed by including the coherent mechanism, the incoherent mechanism and both coherent and 
incoherent mechanisms are indicated by dashed lines, dash-dotted lines and solid lines, respectively. 

• Shift of mean-energy: 

— The incoherent mechanism induces a reduction of the mean energy, while the 
coherent part acts in the opposite way. The origin of this phenomenon can be 
found by looking carefully in figures from |l] to ^. Indeed, we note that both 
coherent and incoherent self-energies induces a shift of the main peaks towards 
lower energy. However, at the same time, a part of strength is shifted toward 
higher energy in the coherent case, which gives rise to a global increase of the 
average mean-energy. Such a behavior can be understood by looking at self- 
energies themselves. An example, in figure |^, the coherent (thin lines) and the 
incoherent (thick lines) self-energies are shown for the GQR in ^°Ca. We see 
that the energy dependence of the real part of the self-energy is different in two 
different mechanisms. While the incoherent mechanism induces a global shift 
of the strength towards lower energies, the real part of the coherent self-energy 
change of sign in the vicinity of the collective energy. This introduces a shift 
toward lower energy of the low energy part of the strength while the high energy 
part is pushed towards higher energies. In some cases, we may even expect that 
a single resonance is splitted into two peaks, as it happens for the GQR in ^°Ca. 

— For the GQR and GMR in lighter nuclei, the shift introduced by the incoherent 
mechanism is stronger that the coherent one. On contrary, the tendency goes the 
opposite direction for heavier nuclei. For the GDR , both effect are comparable. 

— In all cases, the effect of the incoherent mechanism is of the same order of mag- 
nitude as that of the coherent mechanism and can not be neglected in contrast 



to the usual assumption |18 
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• Increase of spreading width: 

— Both the coherent and incoherent self-energy induces an increase of the spreading. 

— In the case of the GQR and GMR, the coherent damping is always much larger 
than the incoherent one while for the GDR both are of same order. For instance, 
in the GQR case, where the strength is in general not Landau fragmented, we 
can see that the coupling to surface modes induces a splitting of the main pic of 
the RPA into different peaks. In the calcium case, this effect can be related to to 
the presence of two collective low-lying 3~ states strongly coupled to the GQR. 

— The fact that the coherent mechanism induces a larger damping can be seen by 
looking at bottom panel of figure 0. In this particular example, we see that the 
imaginary part of the self-energy is larger for the coherent mechanism than for 
the incoherent case, which gives a larger damping width. 

— When the strength is already largely Landau fragmented (like in the GMR of 
^°Ca) both incoherent and coherent effects seems almost negligible. 

— The magnitude of the coherent mechanism becomes larger for heavier nuclei. 

• Additivity of coherent and incoherent effects: 

— The shift in the mean energy, AE, and the increase of the width, AF, are approx- 
imately given by the sum of these quantities obtained by considering the coherent 
and incoherent mechanisms separately. 

The effects of both, coherent as well as incoherent mechanisms appears to be more 
significant for the lighther nuclei as seen from figures 1-4. This may be due to the fact that 
in ligther nuclei a large fraction of nucleons resides in the vicinity of nuclear surface relative 
to the heavier ones, where the predominant effects of both coherent as well as incoherent 
mechanisms occur. Moreover, in light nuclei the energy of the collective states is higher 
leading to an increase of the damping width. 
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5 10 15 20 25 30 

CO (MeV) 

FIG. 7. The real (top) and the imaginary (bottom) part of the coherent (thin hne) and the 
incoherent (thick hne) part of the self-energies for the GQR in ^'^Ca. 

B. Fragmentation of response 

With the high precision experiments, it is possible to determine fragmentation and fine 
structure of the strength functions. In order to characterize systematically the shape of 
strength functions, we present properties of giant resonances spectra in tables |^ - pCV| . 
Depending on the fragmentation of each response, we report average mean energies and 
width for different energy intervals. When the strength is divided into several main peaks, 
we consider energy intervals around the main peaks. Besides the average properties, we also 
report the main peak positions. We emphasize that, in particular for the GQR response, the 
coherent mechanism induces an additional fragmentation of the strength. The incoherent 
mechanism also introduces such a fragmentation, but it is much weaker than the coherent 
effect. The incoherent damping strongly the peak positions. In any case, for a proper 
description of the fragmentation and the fine structure of the strength distributions, both 
the coherent and the incoherent mechanisms should be taken into the description. 
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C. Low-lying states 



The RPA calculations, most often, overestimate the mean-energy of low-lying states. 
We find that the incoherent mechanism reduces systematically the mean-energy of the GQR 
states for medium and heavy nuclei. In table Xll and pCV| , we can see that the mean- 
energy of the low-lying 2+ states is shifted by -1.3 MeV for ^°Zr and ^^°Sn, and by -0.7 MeV 
in ^°^Pb. Such a shift is absent in the calculations with the coherent damping mechanism, 
while it remains when the both mechanisms are included into the description. 



D. Comparison with experiment 

When the strength is highly fragmented, a direct comparison to experimental data is 
hardly possible. Experiments are often analyzed using a fitting plus folding procedure of 
spectra which mix different multipolarities |3^. From this procedure, one extract energy 
(Ea) , width (r^) and percentage of the EWSR {(%EWSR)^). In order to compare with 
experimental data, we convert the experimental data into percentages of the EWSR in given 
energy intervals, which is determined according to 

(25) 

y r^max rQ./27r 

In this expression, a runs over different states of considered multipolarity, and a Lorentzian 
shape is assumed for each mode in the calculations. We note that Gaussian shapes rather 
than Lorentzian only slightly change the reported values. In tables ^ - pCV| , we compare 
the results of our calculations for the percentage of EWSR with experiments in the energy 
interval around the peak energy of the corresponding giant resonance. In all cases, our 
calculations provide a good description for the experiments. In general, the introduction 
of coherent and incoherent mechanisms gives a better description of available experimental 
data. However, in some cases, the percentage of the EWSR obtained in RPA already gives 
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the good order of magnitude. We pay particular attention to the GQR excitations in Ca 
and ^^'^Pb, since both have been extensively studied experimentally and have given long- 
standing discussions 

1. Splitting of the 2+ resonance in '^^Ca 

The GQR response in ^°Ca is known to be split into two components with energy around 
13.5 and 18 MeV with almost an equal fraction of the EWSR (around 30 % to 40 % for both 
peaks). The description of this fragmentation by microscopic calculations is a problem. Only 
recently |[T6| . |17| , microscopic calculations assuming ground state correlations and coupling 
to low lying states reproduce a global splitting. However, these calculations describe the 
global trend of the response and do not provide an explanation for the equal partition of 
the strength. 

Looking at table |VI|, we see that, our calculation with the coherent and incoherent 
mechanism not only reproduce splitting of the strength into two main components (see 
figure [1|) but also give rise to an equal splitting around the main peak (31 % in the interval 
10-16 MeV and 33.6 % in the interval 16-22 MeV), that matches with the experimental 
data. When only the coherent self-energy is included, the calculations can not reproduce 
the splitting, but give the percentage of the EWSR, which are comparable to those obtained 
in ref. i.e. a too high percentage of the EWSR for the second peak and a too low for 
the first one. 

This particular example demonstrates the necessity of taking both coherent and incoher- 
ent damping mechanisms at the same time and illustrates the complementarity of the two 
effects. Indeed, without the coherent mechanism splitting of the strength is not found while 
without the incoherent damping the EWSR is not reproduced. 
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2. Fine structure in ^°Pb 



The calculated strength obtained in our microscopic calculations for the GQR of ^osp]-, 
is displayed in bottom panel of figure |^. As in the case of ^°Ca, global shape of the strength 
exhibits a splitting of the GQR response into two main peaks at 10.8 MeV and 13.8 MeV. 
The first peak is well known and is correctly reproduced by our calculations, while the RPA 
alone does not give a correct description. However the fraction of the EWSR is slightly 
smaller as compared to the experiment (see bottom of table XV ). The second peak has 
never been observed but is also present in second RPA calculations 




9 10 11 
(MeV) 

FIG. 8. Thin line: the strength function for GQR in ^''^Pb, which is calculated including the 
coherent and the incoherent self-energies with a smoothing parameter r] = rjs = 0.025 MeV. Thick 
line: the experimental spectra obtained in the inelastic electron scattering experiment. 



Our calculation, which assume a rather large value of the smoothing parameter only 
gives a global shape of the strength distributions. In order to reveal the fine structure on 
top of the global shape, we also perform calculation s with smaller smoothing parameter 
1] = rjs = 0.025 MeV which corresponds to experimental resolution. The corresponding 
strength distribution is presented in figure ^for the collective energy region 7.6 - 11.8 MeV. 
The calculated response is compared with the inelastic electron scattering data P]. This 
experimental data presents a well-defined fine structure which is also observed with a one- 
to-one correspondence in (p,p') experiments ||70|,|3|]. In figure ^ we see that the calculations 
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agree with the part of spectral properties of the peaks in the vicinity of the collective energy. 
However, we note that below 9 MeV, the fine structure is almost absent in our description. 
The peak positions observed experimentally and obtained in our calculations are reported in 
table pCVl| . We can see from this table that fine structures are already present in the coherent 
case while they are absent in the incoherent one. When both effects are included, it seems 
that part of the peaks are perfectly located as compared to recent (p,p') experiments. It 
has been recently discussed that other peaks might be coming from dipole excitations [|7l] . 
It is also possible that missing peaks might be due to the fact that part of the two-body 
correlations are neglected in the present description or coming from higher order correlations. 



V. CONCLUSION 

In this article, we carry out a systematic investigation the effect of coherent and incoher- 
ent damping mechanisms on the collective response in spherical nuclei at zero temperature. 
Our calculations indicate that both mechanisms play important roles in a proper descrip- 
tion of the nuclear collective response. An extensive comparison with experimental data is 
presented in terms of the fraction of exhausted EWSR for the GMR, GQR and GDR for a 
number of nuclei. We show that the presented calculations are in reasonable agreement with 
the observed collective response. A special attention has been given to the GQR response in 
calcium and lead nuclei where a large amount of experimental and theoretical work exists. 
In particular, we show that, while the usual mean-field theory is unable to explain the equal 
splitting of the 2^ state, the inclusion of both coherent and incoherent damping mechanism 
provides an explanation for fragmentation of the GQR response. Furthermore, by reducing 
smoothing parameter in the calculations, we observe the appearance of fine structure on top 
of the global fragmentation in the strength functions. A comparison with high resolution 
experiment shows that part of the observed peaks energies are located very close to the 
calculated energies. 

Our study demonstrates the importance of coupling to low-lying surface modes for the 
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understanding of fine structures in collective response. For this purpose, the extended mean- 
field description that includes both the incoherent and and the coherent mechanisms in 
an consistent manner appears as a promising tool for the understanding of fine-structure 
properties in the fragmentation of giant resonance excitations. It will be interesting to carry 
our similar investigations at finite temperature. 
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TABLE IV. Properties of GMR in ^°Ca. Top: calculated mean-energy Ei and £"3 and width 
r obtained by integrating moments of the strength in different energy intervals. Calculations are 
carried out within RPA, and by including the coherent, the incoherent and the coherent plus inco- 
herent damping mechanisms, which are indicated in columns under (c), (i) and (c+i), respectively. 
When available experimental data are also reported in the right column. We also display the cal- 
culated width for ijs = WOkeV in parenthesis. Bottom: percentage of the EWSR calculated in 
different energy intervals. In the right column, if available, the corresponding experimental sum 
rules are also reported together with the reactions and references. 
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TABLE V. Same as table |^ for the GDR in ^'^Ca. In addition, in middle panel, Ep^ak indicates 
the positions of the main peaks of the calculated strengths, and the experimental peak position 
and the width of giant resonances are denoted by E and T. 
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TABLE VI. Same as table |V| for the GQR in ^^Ca. 
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TABLE VII. Same as table lIVl for the GMR in 9°Zr. 
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TABLE VIII. Same as table ttVl for the GDR in ^^Zr. 
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23 ± 14 (40Ar+90Zr) || 



TABLE IX. Same as table for the GQR in ^^Zr. 
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i20Sn / 0+ 


RPA 


(c) 


(i) 


Theory 
(c+i) 


Experiment 


Ei-[0-40] 


17.4 


17.8 


17.1 


17.6 




£"3 -[0-40] 


iO.D 


19.7 


18.6 


19.7 




r-[0-40] 


3.6(2.7) 


4.6(4.0) 


3.9(3.1) 


4.8(4.2) 




^peak 


16.8 
18.5 


16.4 


16.1 
17.7 


15.9 


E -15.3 |48] 



3.7 11 



EWSR 



% 



% 



% 



[0 - 40] 
[8 - 20] 

[12 - 20] 



95.9 
90.2 
78.7 



94.7 
65.5 
63.6 



93.8 
77.1 
75.0 



91.7 
65.0 
63.0 



61± 15 (p,p 
72 ± 14 {a, a 
120 {a, a 
74 ± 15 {a, a 
64.6 ± 14 {a, a 



i 

I 

) i 



94. ± 20 (i7O+i20Sn) H 



TABLE X. Same as table Evl for the GMR in ^^OSn 
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i20Sn / 1- 


RPA (c) (i) 


Theory Experiment 

(c+i) 


Ei-[0-40] 
7^0 -\r\ 4(11 

r-[0-40] 


14.8 15.2 14.5 
16 6 18 17 
4.4(3.7) 5.3(4.8) 4.8(4.3) 


14.9 
18.2 
5.6(5.1) 


^-[0-8] 


3.7 3.7 3.7 
5.2 5.3 5.3 


3.7 
5.3 


[10-25] 
E^- [10-25] 


15.2 15.2 14.8 
15.9 16.0 15.6 


14.9 
15.8 


^peak 


2.7 2.7 2.7 
10s 103 103 

±\J.*J ±\J.kJ ±\J.kJ 

14 3 13 4 13 6 
16.9 18.0 


2.7 
10.2 
13.2 

17.9 S=15.4|67] 
r = 4.89 |671 


% EWSR 


% % % 


% % 


[0 - 40] 
[0-8] 
[8 - 20] 
[13 - 18] 


99.0 101.5 97.2 
1.3 1.3 0.8 

85.7 80.6 81.4 

66.1 55.3 57.1 


99.6 
0.8 
78.1 

47.8 62 (7,x) m 


TABLE XL Same as table [V for the GDR in ^^OSn. 
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i20Sn / 2+ 


RPA 


(c) 


(i) 


Theory 


Experiment 










(c+i) 




i?l-[0-40] 


15.8 


16.3 


15.6 


16.1 






1 7 8 


1 Q 6 


1 8 1 


1 q 7 




r-[0-40] 


4.8(4.2) 


6.1(5.7) 


5.3(4.8) 


6.3(6.0) 




Ki - [n si 




5 8 


4 Q 


4 9 




!E^-[0-81 


6.3 


6.3 


5.5 


5.6 




[10-25] 


16.6 


16.6 


16.3 


16.4 




i?3-[10-25] 


17.1 


17.5 


16.9 


17.3 




F, 7 

-'-'peaA: 






4 7 


4 7 








12.8 




12.6 


"K — 1 3 3('0 3^1 Il5^ 












r = 3.7(0.5) m 




16.6 


15.6 


15.7 


15.5 




% EWSR 


% 


% 


% 


% 


% 


fO — 401 


69.8 


68.1 


68.0 


67.0 




[U — 8J 


2 7 


2 7 


2.2 


2 3 




[8 - 20] 


55.1 


43.3 


52.5 


43.0 




[10 - 16] 


40.0 


38.4 


41.5 


38.3 


53 ± 13 (p,p') P 












41 ± 9 (a, a') H 












36 ± 6 (i^O+i20Sn) H 



TABLE XII. Same as table |V| for the GQR in ^^"Sn. 
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208pb / 0+ 


RPA 


(c) 


(0 


Theory 


Experiment 










(c+i) 




i?l-[0-40] 


15.3 


16.1 


15.0 


15.7 




-£•^-[0-401 


16.5 


18.5 


16.8 


18.6 




r-[0-40] 


3.4(2.5) 


4 9^4 4) 


4.0(3.3) 


5.3(4.7) 




-[0-20] 


14.6 


14.4 


14.0 


13.9 






15.0 


14.9 


14.5 


14.5 




r-[0-20] 


2.2(1.4) 


2.4(1.8) 


2.2(1.5) 


2.4(1.9) 






16.3 


14.2 


13.9 


13.6 


E ~13.6 MM 

II 11 ^1 












r ~ 2.5 MM 


% EWSR 


% 


% 


% 


% 


% 




99.2 


102.6 


96.8 


100.1 




[0 — 20] 


88.5 


70.7 


75.9 


65.8 




[10 - 16] 


66.0 


58.4 


68.8 


59.0 


63 ± 17 (d d') 14^ 












64 ± 13 (a, a') |5|] 












61 ± 13 (n n') Il49tl 












59 iHnli 












69 (p,p') m 












64 (3He,3He) |6| 












84 ± 30 (i7O+208pb) m 












92 ± 8 (i7o+208pb) m 



TABLE XIII. Same as table ffvl for the GMR in ^ospb 
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208pb / 1- 


RPA 






Theory 


Experiment 


-Tn 4(11 


13.2 


13.7 


12.8 


13.3 




-C/3-[U-4UJ 


iO.O 


iD.O 


io.y 






r-[0-40] 


4.2(3.5) 


5.2(4.7) 


4.9(4.3) 


5.6(5.1) 




^peak 


2.1 


2.1 


2.1 


2.1 






13.0 


12.5 


11.6 


11.4 






16.7 


16.7 


16.7 


16.7 


E = 13.5 |671 



4.0 |6| 



EWSR 



% 



% 



% 



% 



[0 - 40] 
[0-5] 
[5 - 15] 

[15 - 25] 
[10 - 17] 



99.2 
0.1 
70.7 
23.9 
75.5 



101.8 
0.1 
59.8 
32.3 
66.0 



95.9 
0.1 
65.0 
23.2 
66.1 



97.2 
0.1 
59.2 
28.0 
61.0 



68 (7,x) i 
91 (7,x) i 
90 (7,x) i 
68 (p,p') 13 



TABLE XIV. Same as table HVl for the GDR in ^ospb 
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208pb / 2+ 


RPA 


(c) 


(0 


Theory 


Experiment 










(c+i) 




i?l-[0-40] 


11.5 


12.4 


11.2 


12.0 




-£•^-[0-401 


14.6 


17.6 


15.4 


18.0 




r-[o-4o] 


4 4") 


6.6(6.4) 


5 6f5 2"! 


7.0(6.8) 




Ei-[8-2Q] 


13.2 


13.0 


12.5 


12.6 




i?3-[8-20l 


13.5 


13.9 


13.0 


13.6 




r-[8-20] 




2 8^2 6"! 


1.9(1.5) 


2.8(2.7) 




!Ei-[0-8] 


5.1 


5.1 


4.4 


4.4 






5.9 


5.9 


5.1 


5.1 




r-[o-8] 


1.8(1.6) 


1.8(1.6) 


1.6(1.2) 


1.6(1.2) 




F, 7 

'-'peak 








4 8 






13.0 


11.2 


12.1 


10.8 


'E —10 9 Il53l 












r = 2.4(0.4) 






13.8 




13.8 




% EWSR 


% 


% 


% 


% 


% 


[0 - 401 


92.4 


93.4 


89.0 


90.6 




[0 - 81 


10.7 


10.5 


9.2 


9.1 




[8 — 12.5] 


15.5 


25.2 


39.0 


27.5 


49 ± 12 (p,p') El 












44 ± 10 (p,p') fel 












44 ± 8 I'd n'^l Il69tl 












62 ± 11 (v, d') leiti 












50 ± 8 (0+Pb,H-I) mi 












50 ± 5 (a, a') Q 












33 ± 8 (i7o+208pb) m 












36 ± 9 (20Ne+208pb) H 
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36 ± 12 (40Ar+^u«Pb) p 



208 T 



TABLE XV. Same as table W\ for the GQR in ^ospb. 



(e,e') i (P,P') m (P,P') i (c+i) {c^ 0) 

8.9 8.9 8.9 

9.4 9.3 9.4 9.3 9.3 9.3 

9.6 9.6 9.9 9.9 

10.1 10.1 10.2 10.3 

10.7 10.6 10.7 10.7 10.8 

11.5 11.0 11.0 11.3 

11.9 



TABLE XVL Experimental energies of peaks position observed in the GQR of ^osp]-, (p,p') 
and (e,e') experiments. Calculated peak positions including the coherent, the incoherent and the 
coherent plus incoherent self-energies are reported under columns indicated by (c), (i) and (c+i), 
respectively. 
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